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Pulse responses have been used to calculate frequency response 
functions for several years. Questions have been raised whether use- 
ful pulse responses are obtainable from systems having backlash and 
saturation effects. It has been found in this investigation that by proper 
choice of pulse shape, duration and strength, it is possible to realize 
useful pulse responses despite system defects. 

First it is necessary that the pulse function have continuous deri- 
vatives to avoid infinite inputs. Secondly, the pulse strength must be 
adequate to proyide an output large compared to the backlash and other 
uncertainties. Finally the pulse should endure long enough to introduce 
a pulse of adequate strength without inducing saturation in the system. 

Pulse tests of two simple positional servo - mechanisms were 
performed. The AVSD system is characterized by angular velocity 
signal damping and the MAVSD system by modified angular velocity 
signal damping. Parameters pertinent to the investigation are as follows: 

(1) Backlash 0.25 to 1.0 milliradians 

(2) Gear Train 100 to 1 

(3) Motor FD 84-21-1 Diehl Manufacturing Company 
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(a) maximum acceleration 

(b) maximum velocity 



9300 rad/sec 
415 rad/sec 



2 



3 



30550 



(4) Pulse functions 

(a) Rectangular 

(b) Displaced cosine 

(5) Pulse strength 0.8 to 3.2 milliradian second 

(6) Pulse duration (PD)(NP)R from 0.1 to 1.0 

Rectangular pulse functions were found to be unsuitable because of the 
discontinuous nature of their derivatives. A displaced cosine pulse 
function of (PD)(NP)R = 1.0 and a pulse strength of three point two 
(3.2) milliradian seconds produced negligible response deviations for 
a backlash of one milliradian in both systems. 

Thesis Supervisor: Sidney Lees 

Title: Assistant Professor of 

Aeronautical Engineering 
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OBJECT 

The object of this investigation is to dete rmine a rational basis 
for the choice of pulse functions to test systems with limited non-lin- 
earities. 
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CHAPTER 1 



INTRODUCTION 



Testing individual components and complete systems is a neces- 
sary phase of engineering. This practice not only enables verification 
of predicted performance, but also yields new information concerning 
properties of the system. Prediction of system performance relies on 
the existing knowledge of component and system properties. However, 
new inherent properties are frequently disclosed by the correlation of 
predicted and actual performance. Thus the complementary nature of 
these two aspects of testing is evident. 

Dynamic characteristics of systems and components are inves- 
tigated by analysis of their responses to specific input quantities. Of 
the many possible input quantities, two classes are most useful. These 
consist of transient inducing inputs and sinusoidal inputs. The relation- 
ship between the transient response to a pulse input and the steady- 
state sinusoidal response has been known for many years for systems 
associated with linear differential equations with constant coefficients. 
Practical methods for calculating the sinusoidal frequency response fropn 
the pulse function transient response were devised by Bromberg^ and 
others. 

Pulse function testing is particularly valuable when sinusoidal me- 
thods are impractical. For example, when the life of the system is short, 
as in guided missiles, sinusoidal responses cannot be measured di- 
rectly. Aircraft systems, being subject to changes in environmental 
conditions during test, do not submit readily to sinusoidal response 
measurement. However, the time required for pulse testing is much 
less than that required for sinusoidal techniques. Thus, pulse tests can 
be completed without materially affecting system life and before sig- 
nificant changes in environment occur. In other instances, damage to 
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equipment which can result from sinusoidal tests may be avoided with 
pulse function testing. Recommendations have been made in the interest 
of economy that pulse techniques be practiced in lieu of sinusoidal for 
even those occasions when sinusoidal tests are feasible. 

Objections have been raised to this recommendation. Statements 
have been made <hat physical systems having backlash in the gear 
train together with velocity and acceleration limi’s of the motor will 
yield such poor pulse function responses that the method is useless. 
Simultaneously, it is asserted that this defect is not as important for 
sinusoidal tests. These objections were described by Mr. Rawley 
McCoy of the Reeves Instrument Company, who suggested the subject 
of this thesis. 

From physical reasoning it would appear that pulse techniques 
should provide results at least as good as those obtained from sinus- 
oidal tests in the presence of the non-linearities outlined above. These 
non-linearities exist in some measure in every physical system. 
Successful application of pulse techniques to physical systems are re- 

ported by Seamans, Blasingame and Clementsorr , Bollay^ , Gardner 
(41 

and Ross' ' and others. When difficulties appear because of defects 
in the gear train or limitations of the motor, it seems reasonable to 
suppose that the difficulties may be overcome by the proper choice of 
the pulse function input. The purpose of this investigation is to furnish 
a rational basis for the selection of the pulse function. 

Physical systems are subject to interferences which prevent the 
controlled-member from performing in the ideal manner. Backlash 
in the gear train is comparable to interferences in this respect. System 
design practices are to reduce the effect of interferences to a tolerable 
level described as the level of uncertainty. Similarly, backlash should 
be minimized, preferably to the level of uncertainty, before the system 
performs as desired. The strength of the pulse input must be great 
enough to make the uncertainty and backlash negligible quantities in the 
response. However, when the maximum value of the physical system 
response exceeds a critical value, saturation in the velocity or acceler- 
ation of the controlled-member occurs. The response magnitude is there- 
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fore limited if saturation effects are to be avoided. A range of response 
magnitudes is established by these considerations, between a lower limit 
derived from the backlash or uncertainty level and an upper limit imposed 
by saturation. The problem is to determine whether or not sufficient 
latitude exists in this range to provide necessary information for an 
analysis of the dynamic properties of a given system. In this connection, 
the influence of pulse shape and pulse duration, as well as pulse strength, 
has been examined. A qualitative study of the distortion of the response 
resulting from velocity and acceleration saturation has been made to 
verify the existence of an upper limit on the response magnitude. 

A reasonable surmise is that the deviation from the ideal response 
arising from backlash will be no greater than that due to the general 
uncertainty level. The average deviation created by the backlash is ex- 
pected to be the smaller of the two. This is a consequence of the inter- 
mittant nature of backlash as opposed to the continuous nature of the level 
of uncertainty in a physical system. A consideration of this situation leads 
to the premise that the effect of backlash on the response will be indis- 
tinguishable from the uncertainty level. 

Pulse functions, in this investigation, have been applied to two 
selected systems to provide data for quantitative evaluation of the problems 
associated with pulse testing physical systems. These two systems 
have typical closed-chain structures. Selection of the systems was in- 
fluenced by the Illative simplicity of their performance equations and 
the attendant numerical calculations involved. However, conclusions 
derived from these cases are believed to be quite general in character. 

The investigation reported in this thesis indicates that the res- 
ponses of physical systems derived from pulse tests will be as satisfac- 
tory as those obtained by sinusoidal methods. A definite basis for se- 
lecting the proper pulse function is developed. A pulse shape such as a 
displaced cosine function is desirable. For this type of pulse function, 
it is found that pulse strength can be increased as the pulse duration is 
increased. Thus it is possible to avoid saturation effects and yet res- 
trict the response deviation due to a limited amount of backlash. Under 

these conditions pulse testing of components and complete systems is 

« 

certainly feasible. 
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CHAPTER 2 



DETERMINATION OF PERFORMANCE OF IDEAL LINEAR SYSTEMS 

An ideal linear system is one composed of a set of interacting 
variables so related that a variation of any one variable causes varia- 
tion of one or more of the other variables, and in which the principle 
of superposition applies. Linear systems may be of the open-chain, 
closed-chain, or multiple loop type. A simple closed-chain linear 
system is illustrated in Figure 2-1. 




Fig. 2-1 Functional diagram of a closed-chain linear system 



The performance of an ideal linear system is defined by Draper, 
(5) 

McKay, and Lees' ' as the effects of the system in establishing func- 

-riz 

tional relationships among the physical quantities present. These 
physical quantities fall into two categories, namely system inputs and 
outputs. The inputs in general are independent of actions of the system, 
while the outputs are dependent both upon the inputs and the operation of 
the system. The performance of the system may then be described by 
the functional relationships between the outputs and the inputs. The 
functional relationships are expressed mathematically for linear systems 
in the form of kinematic and dynamic equations of motion of the system. 
These equations take the form of linear integro-differential equations 
with constant coefficients. The coefficients are determined by the para- 
meters of the given system. 
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Draper et al^ develop the concept of the performance function 
which may be derived from the integro-differential equation. The 
functional relationship between a given input and output may be repre- 
sented under zero initial conditions as shown below. 



q out 




q in 



In this expression q Qu ^is the output quantity, q in is the input quantity, 
[lPF| id is the identified performance function. The performance function 
embodies two components. The first, sensitivity, accounts for magnitude 
variations between inputs and outputs of the system. The second compon- 
ent, called the response function, represents the time history of the 
magnitude variations between the inputs and the outputs. 



The form of the performance function of particular interest is 
the steady-state sinusoidal relating function. This form describes the 
response of an ideal linear system to sinusoidal input variations. Sinus- 
oidal inputs may be represented by the rotating vector complex plane 
method or by a plot of the sine function with w^t as its argument and the 
amplitude of the variation as its coefficient. Both methods are illus- 
trated in Fig. 2-2. In this figure <3^ n j a is the amplitude, w^. is the an- 
gular forcing frequency, and t is time. ‘ 








Fig. 2 t 2 Representation of sinusoidal input functions. 
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The output or response of the system will also be a sinusoidal function. 
The response will differ from the forcing function in amplitude and phase 
angle. Dynamic amplitude ratio (DAR) is the ratio of the response am- 
plitude to the input amplitude. Dynamic response angle (DRA) is the 
phase shift between the input forcing function (IFF) and the output 
response (OR). Dynamic amplitude ratio and dynamic response angle 
are both functions of the forcing frequency. Fig. 2-3 shows a typical 
steady-state sinusoidal response for an ideal linear system. The steady- 
state sinusoidal relating function may then be expressed as: 




- (DAR). 



id FV 



j(DRA) id 






WJ 



The steady-state sinusoidal response may be expressed as 



%ut 




eJV 




Fig. 2-3 Representation of the steady-state response of an ideal linear 

system to a sinusoidal input. 
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The variation of the steady-state relating function with the forcing 
frequency is the frequency response of the system. Nondimensional 
performance expressions are derived by dividing the dynamic ampli- 
tude ratio by an arbitrarily selected reference amplitude ratio. Per- 
formance function derivations and frequency responses for first and 
second order ideal linear systems are contained in Instrument Engineering . 
A typical linear scale nondimensional frequency response is shown in 
Fig. 2-4. 
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Figure 2-4 Linear scale sinusoidal response characteristics. 



Pulse functions may be used as forcing functions for ideal linear 
systems. Two types of pulse forcing functions are considered by the 
authors. 

3-1 Rectangular pulse function (RPF) 

Pulse functions of this type are constant for the duration of the 
pulse. Discontinuities exist at the beginning and at the end of the pulse. 
The pulse is described by three attributes. These are pulse height, 
(PH), pulse duration (PD), and pulse strength(PS). To mathematically 
define the rectangular pulse function, let u(t) represent the pulse as 
a function of time, t. 
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u(t) = 0 ; t < 0 

u(t) = (PH) ; 0 - t - (PD) 

u(t) = 0 ; t>(PD) 

Fig 2-5 shows a representation of a rectangular pulse function 




o PD 



TIME 



Fig. 2-5 Representation of a rectangular pulse function. 



Pulse strength (PS) is defined 



r (PD) 

(PS) = (PH) dt = (PH)(PD) 

J o 

.at 

2-2 Displaced ^cosine pulse function (D CPF) 

Pulse functions of this type vary in magnitude throughout the 
duration of the pulse. They attain a maximum value only once during 
the pulse duration. The pulse is described by three characteristics. 
These are pulse height (PH), pulse duration (PD), and pulse strength 
(PS). To mathematically define a displaced cosine pulse function, let 
u(t) represent the pulse as a function of time, t. 



u(t) = 0 ; t < 0 



u ( t) =I™I 
2 



1 - cos 2 r K — — 



(PD) J 



; 0 - t — (PD) 



u(t) = 0 ; t > (PD) 
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Fig. 2-6 shows a representation of a displaced cosine pulse function. 




o 



P D "t ► 



Fig. 2-6 Representation of a displaced cosine pulse function. 
Pulse strength (PS) is defined: 




1 - cos 2W- 

(PD). 



dt 



= (PH)(PD) 

2 



Instrument Engineering ,^) Bromberg^ Seamans, Blasingame, and 
Clementsonf^and others develop numerical and graphical methods for 
obtaining the frequency response of an ideal linear system from data 
obtained through pulse testing. Essentially the methods consist of apply- 
ing a pulse disturbance to the system and recording the input pulse and 
output response as functions of time. Fig. 2-7 shows illustrative time 
responses of second order (0;0, 1, 2) linear systems for pulse function 
inputs. The required performance function form may be determined from 
the time histories recorded. The procedure is to obtain the Fourier 
transforms of the input and the response. These transforms may be 
shown to be continuous functions of frequency. The ratio of the Fourier 
transform of the response to that of the input is the desired relating 
function. The relating function thus determined expresses the complete 
frequency response of the system. 
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Fig. 2-7 Linear second order system time responses to pulse function 

inputs 

From the foregoing it is evident that for an ideal linear system 
the same criteria may be determined by pulse disturbance methods 
as can be determined by sinusoidal methods. This is noteworthy since 
the time required to apply pulse test techniques is much less than that 
required using sinusoidal methods. Further it is conceivable that pulse 
techniques can be applied to systems when sinusoidal testing might cause 
excessive wear or damage to the system whose performance is to be 
determined. 
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CHAPTER 3 



DESCRIPTION OF THE PERFORMANCE OF SYSTEMS WITH LIMITED 

NON-LINEARITIES 



A non-linearity exists in a system when a dependent variable is 
determined by an independent variable in accordance with some non- 
ideal relationship. For example, the independent variable might vary 
sinusoidally and the dependent variable will then vary periodically but 
not sinusoidally. This constitutes non-linear dependence. Fig. 3-1 
illustrates the non-linear relationship just described. 




Fig. 3-1 An example of a non-linear response. 

The non-linearities considered in this thesis are nominal amounts of 
backlash in a gear train as well as saturation characteristics of a 
motor. As noted in Chapter 1, these will exist in some measure in 
almost every physical system. 

A system response forced by a pulse function input is termed a 
transient response. An example of such a response is shown in Fig. 2-7 
of Chapter 2. The transient nature is readily evident. Correspondence 
of the input and output quantities is attained only when both have returned 
to zero displacement. The response does not reach a state of continuous 
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correspondence with the input pulse function prior to this time. Further, 
it does not follow the pulse function with a constant error. The discussion 
which follows is concerned only with system transient characteristics. 

Fig. 3-2 illustrates schematically the relationship existing between 
the controlled-member and the torque summing member when backlash 
is present. The term controlled-member refers to the output member 
of a feedback control system. The torque summing member is considered 
to be the next to last gear in the train. 

Tustin^ develops a technique that makes it possible to closely 
approximate the controlled-member motion when backlash exists in the 
gear train. It is assumed that some point between the torque summing 
member and the controlled-member moves in the manner predicted by 
ideal linear theory. This assumption is reasonable when the backlash 
constitutes a relatively small part of the total motion of the controlled- 
member. The location of the point is established by Tustin's method 
as the mean value of the motions of the controlled-member and the torque 
summing member Motion of this point is referred to as mean motion. The 
actual motion of the controlled-member is calculated by combining the 
mean motion predicted by ideal linear theory with a deviation calculated 
from the backlash. The angular displacement of the point performing 
the mean motion is given by the expression: 

a _ Vm) A (cm) + Vsm) A (tsm) /0 1X 

(mean) 7 (3_1) 

(cm) 1 (tsm) 

where 

= controlled-member inertia 

= torque summing member inertia referred 
to the control member. 



1/ X 

(cm) 

*(tsm) 
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(^P)(tsm) 



TORQUE SUMMING MEMBER 
REFERENCE POSITION 




POSITION 



Fig. 3-2 Schematic diagram of the relationship between controlled - 
member and torque summing member when backlash exists 
in the pear train. 

i 
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The motions of the controlled-member and torque summing mem-* 
ber are intermittant relative to the mean motion as the backlash is taken 
up in one direction or the other. When the backlash is small the motion 
of each member may be represented as a deviation from the mean motion. 
From the definition of the mean motion, it is noted that the deviation of 
the torque summing member is opposite to that of the controlled-member 
when the backlash is being taken up. The sum of these deviations equals 
the whole backlash. Each member may be considered to take up a fixed 
proportion of the total relative motion due to backlash. The portion of 
the backlash taken up by the controlled-member is: 



fysm) A (bl.) 

tbl)(cm) 

l (tsm) T cm ) 



(3-2) 



That taken up by the torque summing member is: 



(cm) A (bl) = A 

~ — A (bl)(tsm) (3-3) 

1 (tsm) (cm) 

It is now possible to describe the motion of the controlled member. 
Referring to Fig. 3 -3a, the decaying transient oscillatory curve repre- 
sents the velocity of the mean point. The backlash is assumed to be 
taken up at the initial instant. As a result, the controlled-member vel- 
ocity duplicates that of the mean point until the first peak occurs. At 
this instant the mean velocity decreases. However, in the absence of 
friction, the controlled-member continues to move with a constant velocity 
because of the backlash. The mean velocity decreases until the effective 
backlash is taken up at which instant the gear teeth make contact. It 
is assumed that contact between the gear teeth is achieved without bounce 
effects. At contact the torque summing member slows down the controlled- 
member while it in turn speeds up. Fig. 3-3b shows the velocity devia- 
tions of the controlled-member relative to the velocity of the mean point. 



30 



POINT OF SEPARATION 




Fig. 3 -3a Velocity of the controlled-member relative to the velocity 
of the mean point for backlash in the gear train. 




Fig. 3 -3b. Deviation velocity of the controlled-member 




Fig. 3-3c. Displacement of controlled-member relative to the 
displacement of the mean point. 



Figure 3-3. 



Motion, of the controlled-member relative to the mean point 
for a decaying transient osciUatiop,of the mean point as 
adapted from Tusti n- '• 
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The controlled-member displacement curve is obtained by integrating 
the velocity curve of the controlled-member. Fig. 3-3c represents 
the angular displacement curves corresponding to the velocities in Fig 3 -3a. 
The integral of each of the shaded triangular areas of Fig. 3 -3a is equal to 
the effective backlash experienced by the controlled-member. 

From the above discussion, it is apparent that the deviation of the 
controlled-member displacement never exceeds the effective backlash 
of the controlled-member. When the backlash itself is small and the 
peak displacement of the mean point is comparably large the net effect 
of backlash is insignificant. For this case, the area between the con- 
trolled-member and mean point angular displacement curves will be a 
very small percentage of the area under the mean point curve. 

Torque saturation characteristics of the motor may be considered to 
impose limits on the velocity and acceleration of the controlled-member. 
The undesirable effects of driving a system into saturation may be seen 
in Fig. 3-4. This figure shows the ideal linear angular displacement 
response of a second order system to a rectangular pulse input. Res- 
ponse curves are also included for both velocity and acceleration 
limiting of the controlled-member. It is noted that velocity and accelera- 
tion limiting increase the peak amplitudes of the transient response. 
Further, a considerable time lag resulting from the limiting is evident. 

Fig. 3-5 illustrates the effect on the velocity of the controlled-member 
of velocity limiting and acceleration limiting. Similarly, Fig. 3-6 
shows the effect of acceleration limiting on the acceleration of the con- 
trolled-member. Certain discontinuities are pointed out in the curves 
of Figs. 3-5 and 3-6. These are attributed to the effect of the discon- 
tinuity in the rectangular pulse function which occurs when it returns to 
zero. These curves were obtained from the General Purpose Simulator 
of the Instrumentation Laboratory of M.I.T. 

Qualitative study of Fig. 3-4 indicates that distortion of the transient 
response resulting from either velocity or acceleration saturation is 
unacceptable. The frequency response calculated by Bromberg's ^ 
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Fig. 3-5 Angular velocity response Of a second order (.o\o t \,z) system 

for a rectangular pulse function input illustrahna +-h© effect- 
of velocity ana acceleration limiting. 
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method from such a transient response will also be distorted. Satis- 
factory correlation between the predicted and actual system response 
under these conditions would be difficult, if possible. It is concluded 
therefore that the pulse strength must be limited to avoid saturation 
effects for successful application of pulse test techniques. 

The character of the pulse shape, as well as the pulse strength, 
becomes important if the controlled- member velocity and acceleration 
are to be prevented from reaching saturation In this connection. Deri- 
vation Summary 3-1 develops the initial accelerations required of repre- 
sentative second and third order systems for a step input. It is evident 
that a step input to the second order system will demand an acceleration 
of the controlled-member at the initial instant equal in magnitude to 
the square of the undamped natural frequency of the system multiplied 
by the magnitude of the step. For the third order system, the same input 
demands an acceleration approximately equal to one and three tenths 
times as great. The argument may be extended to show that a rectangular 
pulse will force such accelerations at both the beginning and end of the 
pulse. Consideration of these facts suggests that acceleration, rather 
than velocity limitations are critical for rectangular pulse functions. 

The problem clearly will not be alleviated by increasing the pulse duration. 

A displaced cosine pulse function, in contrast to the rectangular, 
possesses no discontinuities. Therefore it is anticipated that the relative 
importance of the velocity and acceleration limits will depend both on pulse 
strength and pulse duration. 

The discussion presented here is concerned with a method for cal- 
culating the angular displacement response of a system with backlash in 
the gear train. It is suggested that the response deviations from ideal 
will not exceed the backlash provided velocity and acceleration saturation 
are not encountered. Acceleration saturation will be difficult to avoid for 
rectangular pulse inputs. For small amounts of backlash, it is presumed 
to be possible to obtain sufficiently large peak response amplitudes with 
a displaced cosine input while avoiding saturation. If this can be done, the 
non-linearity of the transient response of the physical system may be 
considered negligible. Close correlation of the predicted and actual system 
responses then becomes practicable. 
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Consider the following second order linear differential equation 
with constant coefficients. 



c\ 

('p^+2 $'w p + 'w ) A. . , = CT 'w 

' ~ W n r \ ' lon+1 



n 



(out) 



n 



(in) 



( 1 ) 



where 

'p = non-dimensional derivative operat or 
A(i n ) = independent variable 

I = dependent variable 

= non-dimensional undamped natural angular frequency 

= damping ratio 

= dimensional similarity factor 

The Laplace transformation yields: 

. 2 
<T 'w 

A, + v (A) = A.. ,(X) (2) 

(out) v/ ' 2 2 ( ln ) v ’ 

(X + 2^'w x+ 'w ) 



*w 

n 

S 

<r 



for X = (X+ j 2^,8 

andX = nondimensional Laplace transform variable 
C( = real component of A 
2 Tf/3= imaginary component of A 



Let A . . x 
(in) 



unit step function 



Then, A (in) (X) = 



(3) 



Substituting in (2) 
A <Out)<*> = 



■w 



n 



)\(A +2.f'w n A+ 'w^ 



(4) 



Derivation Summary 3-1 Initial acceleration of second and third order 
systems for a step input. (Page 1 of 3) 
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A " (out)(A> = A(A 2 + 2fw n X + 'w n :! ) 



<r 'w n A 

A,,(OUt)(A) = (^ + 2tw n X + 'w3 



n 



(5) 

( 6 ) 



For simplicity let 'w = 1.0, 



cr = 1.0 

Substituting for 'w andcn and expanding in partial fractions 



A "<o U «) (A) 






-s-ih-P 



(7) 



2jJ i -S 2 (X-Xj) 2jjl -S 2 (X- A*!) 

Where A 5 ^ is the complex conjugate of Aj . Taking the inverse transform, 

- cf+jfi -:f 2 A'i 1 . -of- j Ji -^ 2 Jfi'* 

2jJi -i 2 



(out) 



2j 



/iAf 



(8) 



where 't = non-dimensional running variable. 
At 't = 0 



A" .. = 1.0 = 'w 
(out) n 



(9) 



(a) initial acceleration of a second order (0;0,1,2) system for a step inppt. 
The third order equation similar to (6) is as follows: 

A ( 1 + 't j X ) 



A" t\\ - 

(out) (A 2 +2.fA+i)(i+'t 2 X) 

when 'w = G~ = 1.0 
n 



( 10 ) 



'tj and 't 2 are first order non-dimensional time constants. 



Derivation Summary 3-1 Initial acceleration of second and third 
order systems for a step input. (Page 2 of 3) 
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CHAPTER 4 



GENERAL DESCRIPTION OF TWO SELECTED PHYSICAL SYSTEMS 

In Chapter 3, a method is proposed for determining the perfor- 
mance of systems with some non-linearities. In order to find out how 
practical the techniques are, two physical systems are selected for 
study. Each is a positional servomechanism with two feedback branches. 
The first is one associated with a second order (0; 0,1,2) performance 
equation. The other is a system associated with a higher order perfor- 
mance equation, namely a third order (1, 0;0, 1, 2, 3) type. Both systems 
have as an input angular displacement A^^ . The output of each sys- 
tern is a controlled-member angular displacement 



4-1 Ideal Linear Properties of the AVSD and MAVSD Systems 

The first system is an angular velocity signal damped positional 
servomechanism which will be referred to as the AVSD system. A 
functional diagram of the AVSD system is shown in Figure 4-1. Its 
ideal linear performance equation is a linear differential equation of 
the following derivative operator form: 



The derivation of the above equation and definitions of the coefficients 
may be found in Appendix A. 

The ASVD system is characterized by the means by which damping 
is achieved. Negative unity feedback of the controlled-member angular 
displacement causes generation of a correction signal whenever the dis- 
placement is different from that of the input. The existence of signal 
level elements in the chain of Fig. 4-1 permits reduction of the correc- 
tion signal by a signal proportional to the torque summing member angular 




intfr 



(4-1) 
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ANGULAR VELOCITY SIGNAL DAMPED POSITIONAL SERVOMECHANISM (psMdVSCl)- 
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FIG. 4-1 FUNCTIONAL DIAGRAM OF BA5IC POSITIONAL SERVO - 
MECHANISM WITH ANGULAR VELOCITY SIGNAL FEEDBACK DAMPING. 



velocity. The reduced correction signal produces an actuating torque 
in the torque generator which tends to drive the correction signal to 
zero. Under these conditions, the actuating torque is the difference 
between a torque proportional to the original correction signal and a 
torque proportional to the angular velocity of the torque summing 
member. Thus a damping-like term, introduced at the signal level, 
is included in the actuating torque. This technique exactly simulates 
viscous damping while avoiding high power dissipation. 

As a final comment on the AVSD system, it can be shown that an 
error will exist between the output and input shafts when the angular 
velocity of the controlled-member is constant. 

The second system studied is shown in the functional diagram of 
Fig. 4-2. The diagram differs from that of Fig. 4-1 by the addition 
of a velocity signal modifier in the angular velocity feedback path. 
Therefore this system will be referred to as the MAVSD system 
hereafter. 

The velocity signal modifier consists of a passive network differ- 
entiating circuit. Incorporating the velocity signal modifier gives the 
following derivative operator form of the system ideal linear performance 
equation. 



This equation is equivalent to a third order linear differential equation 
with constant coefficients. The reader is referred to Appendix A for the 
derivation of the performance equation and the definition of the coeffi- 
cients. 

A frequency dependent damping characteristic results from the 
inclusion of the velocity signal modifier in the MAVSD system. The 
output of the angular velocity signal generator is supplied to the differ- 
entiating network of the velocity signal modifier. The resultant output 




S (fb)(av)[A ; M] P [ A (cm) 
(4-2) 
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FIG. 4-2 FUNCTIONAL DIAGRAM OF 8ASIC POSITIONAL SERVOMECHANISM 
WITH MODIFIED ANGULAR VELOCITY SIGNAL FEEDBACK. DAMPING. 



is then combined with the correction signal generated from negative 
unity feedback of the output position. 

The frequency-dependence of the velocity signal modifier is appar- 
ent in the steady-state frequency response form of its performance func- 
tion. 



With the characteristic time equal tothe undamped natural period, it 
is evident that the function approaches zero for small values of the 
forcing frequency. For values of the forcing frequency of the order of 
the undamped natural frequency, the function approaches unity. The 
undamped natural frequency referred to is that of the second order system 
existing before the inclusion of the velocity signal modifier. 

The desirability of frequency-dependent damping is developed by 
Lees^. Briefly, it is shown in the reference cited that the damper charac- 
teristics should not be prominent for low frequencies, but should become 
important for frequencies in the range of the undamped natural frequency. 
This requirement is met by the velocity signal modifier of the MAVSD 
system. 

4-2 Physical Properties of the AVSD and MAVSD Systems 

So far no consideration has been given to the physical properties 
of these two systems. Generalities have been purposely used to ensure wide 
applicability of possible conclusions. It must be realized that equations 
4-1 and 4-2 are not associated solely with the systems shown in Figs 
4-1 and 4-2. A variety of components may be chosen. For example, 
a large range of sensitivities and gear ratios is available. Equations 
4-1 and 4-2 are still applicable to describe the performance of the systems. 
This demonstrates the power of this method of investigation. 
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To obtain quantitative results from which conclusions may be 
drawn, it is necessary to specify numerical values of the parameters 
involved. Appendix A lists the numerical coefficients of the equations 
of motion of the systems selected for study. These coefficients are 
derived in Appendix A from physical properties of the systems. 

The physical systems chosen use a motor to drive a controlled- 
member in response to an input command. To drive the controlled- 
member, a gear train is often needed. The presence of the gear train 
introduces backlash into the system. In addition interference torques 
or forces are present. The effects of these physical factors must be 
considered. 

It is known that backlash in a system causes a steady state as well 
as a transient displacement error. It is apparent that in a well designed 
system this error must be held to a minimum. To describe this error it 
is necessary to establish a representative value of backlash to be found 
in physical systems. Such a value may be determined by considering 
Fig. 3-2. It is seen that the amount of backlash is the gap between ad- 
jacent teeth when gears are meshed. Angular backlash may be expressed 
in terms of this gap divided by the radius of the gear. Current machining 
tolerance may be specified in ten-thousandths of an inch. Therefore 
the gap may be as large as four times the tolerance. Assuming a toler- 
ance of plus or minus one ten-thousandth of an inch and a gear radius 
of two inches, it is seen that a backlash of approximately one half a 
milliradian is representative. Therefore in this study the range of 
backlash is one -fourth to one milliradian. 

Motors have limitations on output torque, angular velocity, and 
angular acceleration. Each of these limitations must be considered. 
Limited torque specifies a lower limit on the gear ratio that may be 
used. If too low a gear ratio is specified, the mechanical advantage 
of the motor will be insufficient to move the heavier controlled-member 4 . 

A value of gear ratio of one-hundred to one from the motor to the con- 
trolled-member is chosen as typical. The value of gear ratio selected 
limits the controlled -member velocity and acceleration to one one- 
hundredth of the motor maximum velocity and acceleration if saturation 
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is to be avoided. Chapter 3 explains why saturation should be avoided. 

The direct current servo motor with integral tachometer generator 
FD 84-21-1 built by the Diehl Manufacturing Company is considered 
representative of the drive motors available. The maximum speed of 
this motor is 415 radians per second and the maximum acceleration 
is 9,300 radians per second per second. For the selected gear ratio, 
the maximum controlled-member velocity is then 4,150 milliradians 
per second and the maximum acceleration 93,000 miliradians per 
second per second. 

Limitations on controlled-member velocity and acceleration 
establish a maximum error signal the systems investigated may receive 
as an input. In actual use of the systems no problem is presented. This 
is true since the system operates to keep the error signal nulled. In 
pulse testing, however, a limit is established on the pulse strength 
that may be used. Since pulse strength is a function both of pulse height 
and duration, it is possible to introduce too great a pulse height for 
the duration of the pulse. A pulse of this character will cause saturation 
of the drive motor velocity or acceleration or both. 

The relation between the motor inertia and the controlled-member 
inertia is an important parameter. In Chapter 3 it was shown that the 
controlled-member effective backlash is directly proportional to the 
effective inertia oT-the motor and inversely proportional to the sum of 
the effective masses of the motor and the controlled-member. With a 
gear ratio of one hundredth to one the effective mass of the motor when 
considered at the controlled-member is ten thousand times its actual 
value. In this study the effective mass of the motor is made equal to the 
effective mass of the controlled-member. This makes the effective 
backlash to be taken up by the controlled- member one half the backlash 
actually present in the system. Also by this stipulation the actual inertia 
of the controlled-member is ten thousand times the actual inertia of the 
motor. It is evident that this inertia ratio is extreme but it is considered 
desirable. If the actual inertia ratio controlled-member to motor inertia 
is made smaller, the effective controlled-member backlash will also 
be reduced. To emphasize the effects of backlash the ratio of effective 
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inertias for this investigation will be taken as one half. 

Deviations in measurement of physical quantities will result. 

This class of deviations may be described as uncertainties because 
of their unpredictable nature. A maximum range of the output quantity 
may be specified which will include all the uncertainties. This range 
is defined as the uncertainty level. In the systems investigated, the 
uncertainty level is specified to be one half a milliradian. 
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CHAPTER 5 



PHYSICAL SYSTEM RESPONSES 



Physical feedback control systems using a motor and gear train to 
drive a controlled-member must exhibit non-linear performance char- 
acteristics for some operating conditions. For example, when the in- 
put quantities are large enough, the motor will exhibit either acceler- 
ation or velocity saturation or both. Backlash in the gear train will 
cause intermittant non-linear operation when the motion is reversed. 
Chapter 3 outlines a procedure for determining time responses of sys- 
tems with a little backlash. This method will now be applied to the 
physical systems described in Chapter 4. 

The investigation consists of three phases. The first is the cal- 
culation of ideal linear time responses for rectangular and displaced 
cosine pulse input functions. The second part of the study is the deter- 
mination of actdal responses from computed ideal responses when 
backlash is varied from one-fourth to one milliradian. The last phase 
is to evaluate the’ average error resulting from backlash in time re- 
sponses as well as frequency responses of physical systems. 

5rl Assumptions 

The following assumptions are made to facilitate investigation. 

1. The average uncertainty due to interfering forces and tor- 
ques is considered to be zero. As will be shown later this 
uncertainty may be larger than the average error caused by 
backlash. The inclusion of the uncertainties due to interfer- 
ences in this situation obscures the error caused by backlash. 

2. The effects of friction are neglected. The controlled- 
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member velocity therefore remains constant during take-up 

(7) 

of backlash. It is to be noted that Tustin' states that the ef- 
fects of friction are ameliorating. 

3. Bounce effects between driving and driven gears are not 
considered. 

4. Motor acceleration and velocity saturation will be pre- 
vented. To accomplish this, forcing function pulse strenghts 
will be limited. 

5. All initial conditions are zero. The motor gear is posi- 
tioned prior to test in contact with the controlled-member 
gear in the direction the pulse will drive the system when it is 
first applied. Controlled-member and mean displacement re- 
sponses will then be identical until separation first occurs. 

6. The amplitude of the controlled-member displacement is 
considered to be essentially zero when the mean motion dis- 
placement amplitude has settled to zero. This is a necessary 
assumption to establish convergence and integrability of the 
displacement response of the controlle d member. It is a rea- 
sonable assumption when the uncertainty level is comparable 
to the controlled-member effective backlash. 

5r2 Calculation of Ideal Linear Responses 

The Reeves Electronic Analogue Computer (REAC) of the Instru- 
mentation Laboratory, Massachusetts Institute of Technology, was 
used to calculate the ideal linear time responses. The instrumenta- 
tion of the REAC is described in Appendix B. Time responses were 
recorded on the output tables associated with the REAC. Output tables 
were used to obtain the accuracy required for graphical computation. 

Rectangular and displaced cosine pulse functions were applied to 
the AYSD and MAUSD systems. Pulse duration was one-tenth the nat- 
ural period of the AVSD system. Time responses were obtained of 
angular displacement, velocity and acceleration for various damping 
ratios. Angular velocity and acceleration responses were recorded 
for use in determining maximum pulse strengths for each system. 

These responses are shown in nondimensional form in Fig. 5-1 through 
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TIME - undamped natural period ratio T(MP)R = — 
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Controlled -member acceleration ratio 




Time - undamped natural period ratio T(NP)R s jr' 

Fiq:5-9 t>\sPlaced coiin'e pulse fu<\c.t\or\ iripu-t-^nguVar 
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Table 5-1 Numerical values of first peaks of ideal linear 
response of AVSD and MAVSD systems to rectangular 
and displaced cosine pulse functions with 
(PD)(NP)R = 0. 1 and PH = lmr. 
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Fig. 5-9. Numerical values from which these figures were constructed 
are tabulated in Appendix C. Table 5-1 shows the numerical values of 
the response peaks for both systems. 

It may be seen from Table 5-1 that neither pulse function of this 
pulse duration is satisfactory as a forcing function. If a pulse strength 
sufficient to produce a ratio of one hundred to one peak displacement to 
effective backlash is used the resulting motor acceleration will exceed 
its limit and saturation will occur. It is therefore necessary to in- 
crease peak displacement response and investigate the effect on accel- 
eration and velocity peaks. 

Before discussing this part of the investigation, it is considered 

necessary to comment on the acceleration peak values caused by the 

rectangular pulse function. From Derivation Summary 3-1 it may be 

seen that the acceleration peak values should theoretically be indepen- 

2 2 

dent of damping ratio and equal to w^ A jn . w^ for the AVSD system 

is 4007 radians per second per second. It may be seen from Table 5-1 

that the acceleration is dependent on damping ratio for the AVSD system and 

the value is not constant. Further it may be seen from Derivation 

Summary 3-1 that the acceleration peak for the MAVSD system should 

2 

be dependent on damping ratio and should vary from 1. 1 w A to 

2 2 n 

1,4 w A. as damping is decreased from critical, w for the MAVSD 
n in ^ & 2 n 

system is 2800 rad/sec . Referring again to Table 5-1, it may be 
seen that the MAVSD acceleration peaks for damping ratios of 1 . 0 and 
0. 7 are approximately correct. The peak value for a damping ratio 
of 0. 1 is not to be .trusted at all. The cause of this bad reading is at- 
tributed to a scale factor omitted in the computer set-up or subse- 
quent calculations. It must also be pointed out that output table dy- 
namics could have contributed to these errors in computation since 
the acceleration peaks occur at time zero. 

From Derivation Summary 3-1 it may be seen that pulse dura- 
tion has no effect on the acceleration peak caused by a rectangular 
pulse function. This peak always occurs at the instant the system is 
disturbed. Its magnitude is directly proportional to the pulse height. 

It is evident therefore that no decrease in acceleration could result frorq 
increasing pulse duration if pulse height were kept constant. It may be 
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seen from Figs. 5-6 and 5-9 that the maximum acceleration occurs 
later in the response if a displaced cosine pulse is used. It was there- 
fore postulated that increased pulse duration with constant pulse height 
could produce decreased maximum acceleration if the forcing function 
were a displaced cosine pulse. 

The investigation was then continued using only the displaced co- 
sine pulse function as the forcing function. Only the seven-tenths 
critically damped form of each system is studied. Pulse duration is 
varied from one -tenth to one natural period of the AVSD system. The 
calculated responses in non-dimensional form are shown in Figs. 5-10 
through 5-15 for both systems. Appendix C tabulates the numerical 
values from which these figures were constructed. From these fig- 
ures it may be seen that for constant pulse height as the pulse dura- 
tion increases 

1. Angular displacement peak increases. 

2. The first angular velocity peak increases until the pulse 
duration is one half the natural period. Thereafter it 
decreases. 

3. The second angular velocity peak increases. 

4. The first angular acceleration peak decreases. 

5. The second angular acceleration peak exceeds the first 
angular acceleration peak when the pulse duration is 
greater than three -tenths the natural period of the AVSD 
system. 

6. The second angular acceleration peak decreases after 
the pulse duration exceeds three-tenths the natural 
period. 

The investigation so far has shown that it is possible to increase ti 
the controlled member angular displacement peak while decreasing the 
acceleration peak. 

Table 5-2 lists the numerical values of the response peaks of the 
AVSD system to a displaced cosine pulse function for various values of 
pulse duration. Table 5-3 lists the same information for the MAVSD 
system. From these tables it may be seen that the desired ratio of 
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Time - undamped natural period, ratio T(HP)R, 

Fig. 5-12. Displaced cosine pulse function in put- angular acceleration response 
curves for trie AVSD system for variation of the pulse duration. (.PwJ) 
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Table 5-2 Numerical values of peaks of ideal linear responses 
of AVSD system DR = 0.7 to a displaced 
cosine pulse function. 



PS 


PD 


A 

cm 

mr. 


• 

A - mr/ 
cm ' sec 


A -mr/ 2 
cm ' sec 


PS/ 

(ref) 


NP 


1 st Peak 


2 nd Peak 


1 st Peak 


2 nd Peak 


1.0 


0 . 1 


. 1449 


14. 30 


1 

4.31 


3305 


1502 


3. 0 


0. 3 


.419 


25. 35 . 


12. 12 


2277 


2480 


5. 0 


0.5 


. 633 


27. 24 


18. 34 


1641 


2384 


7. 0 


0.7 


. 771 


26. 15 


22. 04 


1249 


2004 


10. 0 


1.0 


. 891 


23. 54 


23.25 


858. 5 


1426 


1.0 


0 . 1 


. 1449 


14. 30 


4.31 


3305 


1502 


1.0 


0. 3 


. 1397 


8.45 


4.04 


759 


826. 6 


1.0 


0. 5 


. 1266 


5.45 


3. 67 


328.2 


476. 8 


1.0 


0.7 


. 1101 


3. 74 


3.15 


178.4 


286. 3 


1.0 


1.0 


.0891 


2.35 


2.33 


85. 85 


142. 6 



PS 



(ref) 



0. 00993 mr sec. 
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Table 5-3 Numerical values of peaks of ideal linear responses 
of the MAVSD system DR = 0. 7 to a displaced 
cosine pulse function 



PS 


PD 


A 

cm 

mr. 


A - mr/ 
cm ' sec 


A -mr/ 2 

cm ' sec 


PS/ 

(ref) 


NP 


1 st Peak 


2 nd Peak 


1 st Peak 


2 nd Peak 


1 . 0 


0 . 1 


. 1485 


14. 12 


3. 75 


3215 


1343 


3. 0 


0. 3 


.432 


25. 25 


11. 01 


2250 


2300 


5. 0 


0. 5 


. 657 


27. 8 


16. 85 


1679 


2282 


7.0 


0. 7 


.818 


27. 5 


21. 0 


1257 


1930 


10.0 


1.0 


.967 


25. 0 


23. 9 


888 


1441 


1.0 


0 . 1 


. 1485 


14. 12 


3. 75 


3215 


1343 


1.0 


0. 3 


. 1472 


8.43 


3. 67 


750 


767 


1.0 


0. 5 


. 1315 


5.56 


3. 37 


336 


457 


1.0 


0. 7 


.1169 


3. 93 


3. 00 


179. 6 


276 


1.0 


1.0 


.0967 


2.50 


2. 39 


88.8 


144. 1 



PS(ref) = 0. 00993 mr sec. 
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angular displacement peak to effective backlash may be achieved. 

Pulse durations of seven-tenths the natural period of the AVSD system 
or greater may be used. Smaller ratios of displacement peak to ef- 
fective backlash may be obtained by using shorter pulse durations or 
smaller pulse strengths at the large values of pulse duration. 

5-3 Calculation of Actual Responses 

The second part of the investigation deals with the extrapolation 
to the actual response curves of the AVSD and MAVSD systems. This 
is accomplished using the method outlined in Chapter 3. The major 
portion of this phase is graphical integration and construction. In 
Appendix D will be found the complete calculation of a typical control- 
led-member angular displacement response curve. This section will 
consist only of a detailed description of the application of the method 
outlined in Chapter 3. 

The starting point is the ideal linear response peak table. Know- 
ing the acceleration and velocity limits of the motor and the gear ratio 
of the physical system to be tested, these limits are reflected through 
the gear ratio. These limits then become the maximum controlled - 
member velocity and acceleration permissable. These values for the 
AVSD and MAVSD system are respectively 4, 150 milliradians per sec- 
ond and 93, 000 milliradians per second per second. The next step is 
the determination of the maximum pulse height that may be used to 
disturb the system. In this investigation it can be seen that the accel- 
eration limit is the deciding limit. All calculations are made for a 
pulse duration equal to the natural period of the AVSD system. 

To determine the maximum pulse height, divide the allowable 
controlled-member velocity and accelerations by those tabulated in 
Tables 5-2 and 5-3 for the selected pulse duration. It may be seen that 
this maximum pulse height is 65.2 milliradians for the AVSD system 
and 64. 5 milliradians for the MAVSD system. The mean member ang- 
ular displacement response may now be computed. This is accom- 
plished by dimensionalizing Fig. 5-10 for the AVSD system and Fig. 
5-13 for the MAVSD system is such a manner that the peak displace- 
ments are respectively 58.2 and 63.1 milliradians. The solid curve 
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in Fig. 5-16 represents the mean member motion of the MAVSD sys- 
tem. The next step is to dimensionalize Fig. 5-11 for the AVSD sys- 
tem and Fig. 5-14 for the MAVSD system. This is accomplished by 
making the mean member velocity first peaks 1524 milliradians per 
second and 1612 milliradians per second. The solid curve in Fig. 5-17 
represents the mean member velocity of the MAVSD system. The 
next step is to determine the controlled-member velocity curve. It 
will be recalled that friction effects were neglected in this investiga- 
tion. Therefore the controlled-member velocity is identical with 
mean member velocity until separation occurs. Separation occurs 
when the mean member velocity is a maximum or a minimum. At 
this instant of time the controlled-member velocity is at the maximum 
or minimum mean member velocity. The controlled-member contin- 
ues at this velocity until the area between the mean motion velocity 
curve and the controlled-member velocity curve equals the controlled- 
member effective backlash. At this instant the controlled-member 
velocity assumes the mean member velocity. The controlled-member 
velocity curve for the MAVSD system with one quarter milliradian of 
effective backlash is shown in Fig. 5-17 as the solid curve including 
the shaded triangular areas. The hatched areas represent one quarter 
milliradian. The dashed curve of Fig. 5-16 is constructed by plotting 
the area of the hatched area of Fig. 5-17 as the ordinate difference 
from the solid curve of Fig. 5-16 from time of separation until the 
two velocity curves are identical. It is convenient to use a difference 
plot to show the variation of the difference ordinate. The difference 
plot for the MAVSD system is shown in Fig. 5-18 for an effective con- 
trolled-member backlash of one fourth millirad ian.This difference is 
defined as (£)f (t). 

The method of determining the controlled-member angular dis- 
placement and velocity responses for one value of pulse height and 
one value of effective backlash for the MAVSD system has been ex- 
plained and illustrated. These steps must be repeated for each value 
of effective backlash at each pulse strength desired. No curves are 
included for other values of pulse height or effective backlash. 
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Fig. 5-16 Plot of angular displacement of mean member and controlled-member of 
MAVSD system for displaced cosine pulse function input showing effect of backlash. 
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Time - 1 - seconds 

Fig. 5-17 Plot of angular velocity of mean member and controlled member of MAVSD 
system for displaced cosine pulse function input showing effect of backlash. 
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0 0.02 0 . 04 - 0.06 0-08 0.\0 0A7, 

Time - t • seconds 

Fig. 5-18 Plot of controlled-member angular displacement deviation due to backlash 

in the MAVSD system with DR *0.7 for displaced cosine pulse function input 
(PDJNP)R * 1.0 



5r4 Calculation of Average Percentage Error 



The final phase of this investigation consists of the calcualtion of 
the average percentage error in magnitude of the angular displacement 
of the controlled-member time response and frequency response. 
These functions are defined and derived in Derivation Summary 5-1. 



The angular displacement time response of a physical system 
with limited nonlinearities may be written in the form 

ADTTR (pg) = f(t) + C)f(t) (1) 

where ADTTR^ pg j = angular displacement transient time response 

of the physical system 

f (t) = angular displacement transient time response 
of the ideal linear system with same paramet- 
ers as the physical system but with no non- 
linearities 

(£)f(t) = angular displacement deviation of ADTTR^ s j 

from f (t) 

Each quantity defined above is a function of time (t). <£)f (t) is calcu- 

lated as outlined in Chapter 3. 

(a) Angular displacement time response 

— 

Cramer shows that Schwarz inequality may be written in the 
form 



J |(f)t+(£)f(t)| 2 dt<J [|f(t)| 2 +2|(£)f(t)||f(t)|+|(£)f(t)| 2 jdt (2) 



,00 



Neglect the term f | (£)f(t) 

J -00 



dt as small compared to other terms. 



Derivation Summary 5-1. Derivation of error describing functions for 
a physical system with limited non-linearities. (Page 1 of 5.) 
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Subtract 




from both sides of equation (2) 



_°o „ n.00 I oo 

\ |f(t)+(£)f<t)| Z dt- ( |f(t)| Z dt£( 1 2 (£ )f (t)| I f(t)|dt (3) 

J -00 vJ -oo ^ -00 

o 2 

Divide Equation (3) by l | f(t)l dt and multiplying by 100% 

J -oo 



~0O 2 9 ^°° 

100% j |f(t)+(£)f(t)| dt-100% j I f(t)rdt 100% j 2|(£.)f(t)||f(t)|dt 



r i f(t, i 

i ) — 



dt 



-oo 



r°° i 

i |t(t> 



dt 



-00 



(4) 



Define 



APEM 






100 f |f(t)+(£)f(t|| 2 dt -100 f | f (t)| 2 dt 
J -oo J -00 



(ttr) 



(5) 



f M 1 



dt 



where 

APEM^ r j '*= average percentage error in amplitude of the tran- 
% sient time response of a physical system 

it 

Therefore 

r* 00 

200 ( |(£)f(t)||f(t)| dt 

J -00 



APEM,, , . < 
(ttr) — 



/>°° 0 

V f(t) A dt 

J -CO 



in percent (6) 



(b) Expression of average percentage error in amplitude 
of the transient time response of a physical system 



Derivation Summary 5-1 Derivation of error describing functions for 
a physical system with limited non-linearities. (Page 2 of 5.) 
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The angular displacement frequency response of a physical sys- 
tem with limited non-linearities may be written' in the form 



ADFR (pg) = F 0) + (£) F (0) 



(7) 



where 



ADFR^ ps j = angular displacement frequency response of 
the physical system 

F (|3) = angular displacement frequency response of 

the ideal linear system with same parameters 
as the physical system but with no non-linear- 
ities 

(£) F (fi) = deviation of ADFR^ pg ^ from F (|3) 

|3 = frequency ratio and is the forcing frequency 

divided by some arbitrary reference frequency 

(c) Angular displacement frequency response 
Applying the steps outlined in section (b) it can be shown that 



100 J | F(3)+(e>F(3)| 2 d3 - 100 ^ |F(3)| 2 d(3 200j |(£>F(|3)||F</3)|d3 



/>°° 

\ |f (3)| 
J -00 



d£ 



f |f o >| 

J -00 



d3 



( 8 ) 

The mathematical expression F(0) is used to describe a physical re- 
lationship. In nearly all physical systems F(/3) is an even function. 
Negative frequencies do not exist. Therefore 



Derivation Summary 5-1 Derivation of error describing functions for 
a physical system with limited non-linearities. (Page 3 of 5.) 
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100 ^ |F(0)+(£)F(/3)| 2 d^lOO J |F (/3 ) | 2 d0 400 |(£)F0)| |F(|3)|d|3 



(9) 



\ |FO)| 

J - oo 



Define 



APEM 



A 



'dj8 2yjFO)| 2 

100 f |F(/3)+(£)F(/3)| 2 d/3-100 f |F(j3)| 2 d|3 
J -00 J -00 



(fr) 



( 10 ) 



f M : 

J -00 



d/3 



where 

APEM^^j = average percentage error in amplitude of the fre- 
quency response of a physical system 



Therefore 



200 f°°|(£)F0)||F(|B)|d0 

Jo 



APEM (fr j< 



in percent (11) 



f |f(3)| 

Jo ' 



d(3 



(d) Expression of average percentage error in amplitude of 



the frequency response of a physical system 

,2k. 

Morse and Feshback^^ show that Parseval's formula may be 
written 






( 12 ) 



From Equation (12) it may be seen that 



Derivation Summary 5-1 Derivation of error describing functions for 
a physical system with limited non-linearities. (Page 4 of 5. ) 
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f°° |F(j3)+(£)F^)| 2 d|3 - r°° [F(J3)| 2 d0 f |f(t)+(£)f(t>| 2 dt- f |f(t)| 2 dt 
J -00 J -00 J -00 J -oo 




d/3 




dt 



Therefore 

APEM (fr) = APEM (ttr) 

(e) Relationship between average percentage error in amp- 
litude of the transient time response and average per- 
centage error in amplitude of the frequency response of 
a physical system 



Derivation Summary 5-1 Derivation of error describing functions for 
a physical system with limited non-linearities. (Page 5 of 5. ) 



The deviation plots described in 5.3 are plots of (£)f(t) as stated. 
Plots similar to the solid curve of Fig. 5-10 are constructed for each 
value of pulse strength to be considered for each system. A plot of 
2| (f)f(t) ||f(t)|is then constructed for each effective backlash value by 
multiplying the ordinates of the (£)f(t) and f(t) curves together. A plot 
of |f(t)| 2 is constructed for each pulse height for each system. These 
curves are integrated. Graphical methods were used by the authors. 
By combining the integrated functions as shown in section of Deriva- 
tion Summary 5-1 APEM^.^ is obtained. 

Fig. 5-19 is a nondimensional plot of APEM^ tr j as a function of 

pulse strength and effective backlash for the AVSD system. Fig. 5-20 
is a similar plot for the MAVSD system. The quantities used to non- 
dimensionalize these curves are effective backlash ratio-BR- and ref- 
erence pulse strength -( PS) re ^ These two quantities are defined be- 
low. 
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Average percentage error in transient time response APEM(-ttr) 




O O.z 0.4 0.6 0.8 1.0 



Pulse 



Strength. Ratio 



PS 



Fig. 5-19 Angular displacement average error produced by backlash 
in the AVSD system with DR=0. 7 and (PD)(NP)R = 1.0 
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Fig. 5-20 Angular displacement average error produced 
by backlash in the MAVSD system with (DRV* 
0. 7 and (PD)(NP)R s 1.0. 
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BR 



4 kl(cm) n 




T 

n 



where 



bl 



(cm) 
n 




controlled-member effective backlash 

gear ratio motor to controlled- member 
maximum motor velocity 
natural period of AVSD system 



(PS) 



ref 




(PH) 



max 



( 1 



_ 2tt . x 
cos t ) 

n 



dt 



where (PH) = pulse duration 

(PH) maJ£ = maximum pulse height possible avoiding motor 
saturation 

t = time 

= natural angular frequency of AVSD system. 



Table 5-4 lists the numerical values from which Fig. 5=19 was con- 
structed. Table 5-5 tabulates the same information for Fig. 5-20. 
From Derivation Summary 5-1 it may be seen that APEM^ r ^ 

is equal to APEM^.^ . Therefore although not so labelled. Figs. 5-19 

and 5-20 show the variation in APEM /r * for different values of BR 

(fr) 

and (PS) /pg. for the AVSD and MAVSD systems respectively. 

' a 'ref 

The effect of the uncertainty level in a system is handled in an 
analogous manner. It is apparant that the uncertainty level is present 
continuously throughout the motion. For this case, (£)f(t) remains 
constant and is equal to the magnitude of the uncertainty level (U)f(t). 
To compare the effect of uncertainties and backlash, a sample com- 
putation of the average percent error due to an uncertainty level of 
one half milliradian in the MAVSD system is included in Appendix D. 
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Table 5-4 Average percentage error in transient time response 
of the AVSD system due to backlash. 



PSR 

A (bl)(cm) 

(mr)\ 


1.0 


0. 75 


0.5 


0. 25 


APEM 


APEM 


APEM 


APEM 


0. 500 


1. 676 


2.22 


3. 32 


6. 59 


0. 375 


1.25 


1. 67 


2.49 


4. 94 


0. 250 


0. 83 


1. 11 


1. 66 


3. 30 


0. 125 


0.42 


0.55 


0. 83 


1. 65 



Table 5-5 Average percentage error in transient time response 
of the MAVSD system due to backlash. 



X. PSR 


1.0 


0. 75 


0. 5 


0.25 


A (bl)(cm) 

(mr)v^ 


APEM 


APEM 


APEM 


APEM 


0. 500 


1.550 


2. 065 


3. 10 


6. 20 


0. 375 


1.202 


1. 602 


2.404 


4.81 


0. 250 


0. 798 


1. 064 


1.596 


3. 19 


0. 125 


0. 394 


0. 525 


0. 788 


1.576 



84 



The results indicate that an average error of 2.42 percent will exist 
due to the uncertainty level. Fig. 5-20 shows the average error for a 
total backlash of one half milliradian to be 0. 798 percent for a pulse 
strength ratio of 1. 0. It is concluded from this comparison that when 
the magnitudes of the backlash and uncertainties are equal, the re- 
sponse deviations due to backlash are negligible. 

Investigation of the technique of pulse distrubance testing of feed- 
back control systems discolses several significant considerations. 
Velocity and acceleration limitations of a given system must not be 
exceeded. The peak amplitude of the displacement response of the 
system must be large relative to the backlash. These requirements 
can be met simultaneously by proper choice of the pulse function. In 
this connection, pulse shape and duration, as well as pulse strength, 
are important. Utilization of the correct pulse function minimizes 
the response deviations due to backlash, and these deviations are neg- 
ligible for nominal amounts of backlash. 
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CHAPTER 6 



CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDY 

Responses of physical systems to pulse functions are used for 
evaluating the linear properties of the systems. The intent of this 
study is to determine how to apply pulse function inputs which yield 
information about the linear properties of systems having some signi- 
ficant non-linearities. 

Physical systems may exhibit a wide range and variety of non- 
linearities. In order to develop methods for studying such systems, 
the non-linearities have been restricted to backlash in the gear train 
and saturation of the torque source. Velocity saturation exists when the 
torque source output member has a maximum velocity it cannot exceed. 
Similarly, acceleration saturation implies the existence of an upper 
acceleration limit. These are quite general since every physical systerp 
possesses maximum velocity and acceleration limits. The procedures 
described in the previous chapters are adaptable to other types of non- 
linearities. 

Two representative physical systems have been studied here. 

The AVSD system is a second order type, characterized by angular 
velocity signal feedback damping. Modifying the angular velocity feed- 
back signal converts it into the third order, MAVSD system. Non- 
linearities in the form of small amounts of bac klash as well as velocity 
and acceleration limits on the motor are introduced in the ideal systems. 
Rectangular and displaced cosine pulse function inputs are applied and the 
system responses studied. The effects of pulse shape, duration and strength 
on the system responses are examined. The following conclusions are 
derived from this study. 
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6-1 Conclusions 



(a) Velocity and acceleration limits of the motor must not be ex- 
ceeded Saturation of any component can distort the transient response 
beyond usable limits. Correlation with the ideal responses determined 
from linear analysis is then impossible. 

(b) The pulse function input must have continuous derivatives to 
avoid infinite inputs. This is first described in the theoretical analysis 
of Derivation Summary 3-1 of Chapter 3. The components of a physical 
system are incapable of accepting an infinite input. They may be driven 
into saturation by such an input. 

(c) The pulse strength must be adequate to force an output large 
compared to backlash and uncertainties. A ratio of the peak amplitude 
of the transient response to the effective controlled- member backlash 
of one hundred to one is desirable. If this is achieved without saturation 
distortion, the deviations from the ideal response are negligible. 

(d) The pulse duration should be approximately equal to the undamped 
natural period of the dominant mode of the system. Under this condition 
the pulse strength may be adequate without inducing saturation. 

(e) A rectangular pulse function is unsuitable for the systems 
in this study. This is discussed from a theoretical point of view in 
Derivation Summary 3-1 of Chapter3 and verified numerically in Chapter 5. 
The discontinuous nature of this function and its derivatives makes it 
impossible to avoid saturation of the acceleration and velocity for the in- 
vestigated systems. This condition is true regardless of pulse duration. 

(f) Functions such as the displaced cosine pulse may meet the 
requirements stated above. They possess continuous derivatives. Further, 
it is possible to introduce adequate pulse strength by adjusting the pulse 
duration without inducing saturation. 

(g) The average non-dimensional error in amplitude due to backlash 
and uncertainties is the same in both the time and frequency domains. This 
is shown in Derivation Summary 5-1 of Chapter 5. Since the correct 
choice of a pulse function results in a negligible average amplitude error 
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in the transient response, the average error in the amplitude of the 
calculated frequency response will also be negligible. 

(h) When the backlash and the uncertainties are of the same order 
of magnitude, the average amplitude error due to backlash is smaller 
than that due to the uncertainties. 

6-2 Recommendations for further study 

The investigation reported in this thesis is limited to a study of the 
non-dimensional error in the frequency response amplitude. The effect 
of backlash, uncertainties and saturation upon the dynamic response 
angle has not been determined. However, the necessary information 
can be obtained by an extended study of the kind used herein. It is 
recommended that this be done. 

Further study of the problem of calculat ing the frequency response 
from the transient response is suggested. Cognizance of the character 
of transient response errors arising from backlash and other uncertainties 
is important here. The distribution of errors in the frequency response 
may be ascertained from carrying out this recommendation. 

A theoretical investigation of the sinusoidal responses of systems 
having backlash and saturation effects will be useful in comparing direct 
frequency response testing methods with pulse techniques. 

-V 

Experimental verification of the conclusions described above is 
strongly recommended. More convincing conclusions may then be drav/n 
as to the relative merits of pulse techniques versus direct frequency 
response testing methods. 
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APPENDIX A 



DERIVATION OF THE PERFORMANCE EQUATIONS OF THE SYSTEMS 



The performance equations of the positional servomechanisms 
described in Chapter 4 are derived by summing the torques acting on 
the controlled-member. 

A - 1 . Positio nal Servomechanism with Angular V elocity Feedback 
Damping .-A VSD System 

Consider the components of the functional diagram of Fig. 4- 1. No- 
tation used is defined there. 

(a.) (Sc) = angle signal comparator 



INVESTIGATED 





(A - 1) 



(b.) (pa) = Preamplifier 



e (Pa) S (pa) [e;ej e (sc) 



(A -2) 



(c.) (vsc) = Velocity signal comparator 



e 



(vsc) 




(A -3) 



(d.) (da) = Drive amplifier 



e 



(da) 



S 



(da) [e;ej e (vsc) 



(A -4) 
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(e.) 

(tg) = Torque Generator 

M (tg) = S (gt) [A om ; A tsm J S (tg)[e; Mj e (da) 

2 * 

" S <8** [ A cm : A t sm ] S (tg)lA;Ml A (cm) (A-5) 

Reduction of the output of the torque generator indicated in equation (A-5) 
is the result of motion of the torque generator rotor relative to its case. 

(*•) 

(avsg) = Angular Velocity Signal Generator 

e (avsg) S (avsg)£A;e3 A (tsm) 

2 4 

6 (avsg) |[ A C m ; A tsmJ S(avsg) A ( cm ) (A ~ 6 * 

(g) 

(vsa) = Velocity Signal Amplifier 

e (vsa) ^(vsa) [e;ej e (avsg) ( A ~ 7 ) 

Combining the above relationships for the inner and outer loops of the 
functional diagram. Fig. 4 - 1 . yields: 

M (tg) S (gt)[^cm ;A tsm| S (da)[e;e] S (vsc) [e;e] S (pa)[e;e| 

S (sc)[A;e]< C > A (cm) ' S (gt) [A cm ; A tsm ] S (avsg)[A;^| S (vsa)[e;e] 

S (vsc)[e;e] S (da)[e ; e] S (tg) [e;M] A (cm) ' ^gt)^ cm ;A tgm j ^)[a;m] ^(cm) 

(A -8) 
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An inertia reaction torque is present when the controlled-member is 
accelerated with respect to inertial space. 



l (ir) \cm)(eq) A Tl - cmj ^ 

If the base does not move with respect to inertial space. 



A |l - cmj ^(cm) 

M, . . = - I, . A / , 

(lr) (cm)(eg) (cm) 



(A -10) 



A basic dampingiorque exists due to the relative motion between the con- 
trolled-member and the supporting base. It is partially due to friction 
in the bearings and gear train, and to air damping. 




C d A 



0 - 



cm> 



This reduces to 



(A - 1 1) 



M 



(d) 



= - C , A. . 
d (cm) 



(A- 12) 



if the base does not move with respect to inertial space. 



Torque Summing Equation at the controlled-member 



M. . + M,. . + M. + M,. ... =0 (A-13) 

(tg) (lr) (d) (intfr) v ' 

Substituting equation (A-8), (A-9), and (A-12) in equation (A-13) gives 
the following performance equation of the AVSD system; 



(gt)[A cm ;A tsm] S (vsc)[e;e] S (pa) [e;e] b (sc) |A;ej (C) A (qm) 

(gt) |A C m ;A tem| ^ avs S)[ A ' e ] S (vsa)[e;e] S (vsc)re:e] S (da) [e;e] ^tg)[e;M] 

- fcm)(eq) A(Cm) + M (intfr. ) = ° < A ' 14 > 



S 

t 

+s 






+ C. 



A A 

V 



cm 
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Define: 



S ( P s) A;M = S (gt)[A cm ;A ts . m J S (tg)[e; M] S (da)[e;e| S (vsc)[e,eJ S (pa)[e;e] 

S( sc) [*A * ej = e ^ as ^^ c coefficient (A-15) 

S (fb)(av)[A;JVl] " S (avsg)[A;e] S (vsa)[e;e] S (vsc) [e;e| S (da)[^e| S (tg)[e;M| 

(A - 1 6) 

S ( P s)[A;m] = S (gt) 2 [A cm ;A tsm ]^ S (tg)[A;M] + S (fb)(av)[A;]VlJ + C d 

(A - 17) 



Substitution of equations (A -15) through (A -17) into equation (A -14) 
and rewriting the result in operator form yields: 




(cm)(eq) 



P +S (ps)[A;M] 



^(cm) S(ps)£A;Mj^^(cm) + ^(intfr) 

(A- 18) 



A -2. Positional Servomechanism with Modified Angular Veolcity Feed- 
back Damping - MAVSD System . 

Consider the components of the functional diagram of Figure 4-2. 
Note that the only difference in components from those of Figure 4-1 
is the addition of a velocity signal modifier in the angular velocity feed- 
back path. The velocity signal modifier is a simple passive network 
differentiating circuit which has the following performance function: 



PF 



(vsm)[e;ej 



(CT) (vsm) p 
l + < CT *vsm)P 



(A-19) 



With the velocity signal modifier in the angular velocity feedback path 
the torque generator equation is as follows: 
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M (tg) " S (ps)[A:M] (C)A (cm) " 1+iCTl p 

(L, Mvsm) P 





vsm) 



S 



(fb)(av)[A;M] A (cm) S (gt)[A 



A cm ;A tsm3 



1^tg)[A;M] A (cm) 



(A -20) 



The inertia reaction and damping torques derived in equations (A -10) 
and (A -12) are unaffected by the velocity signal modifier. Subsituting 
the expressions for the torques into the torque summing equation (A- 13), 
simplifying and writing in operator form produces the performance 
equation of the MASVD system. 



Attention is invited to the fact that the coefficients of equations (A -18) 
and (A-21) may be kept constant while permitting the sensitivity of the 
gear train to vary, This is accomplished by adjusting amplifier sensi- 
tivities where necessary to compensate for a change in gear train 
sensitivity. As a consequence, it is possible to use the gear train 
sensitivity as a variable parameter in the study of pulse testing techniques 
lor systems with backlash. 

A -3 Numerical Performance Equations for the A VSD and MAVSD Syst ems 

Numerical values for the coefficients of the performance equations 
derived above were selected from Lees^ 10 ^ for a system damping ratio 
of seven-tenths critical at high frequencies. Coefficients were then 
determined for critical damping and one-tenth critical damping without 
changing the undamped natural frequency. Interference torques were not 



I 





^fb)(av)[A;]V^ P ^cm) 



1 + < CT \vsm)P 



= S (ps)[A:Mj (C ) A (cm) + M (intfr) 



(A-21) 
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considered in this investigation. These equations are shown below. 



AVSD System - undamped natural frequency w n ^p S ^ = 63.3 rad/sec 

(1) Critical Damping, (DR)^p g ^ = 1.0 

(p 2 + 126. 6p + 4007) A (cm) = 4007 A ( - n) (A-22) 



(2) Seven-tenths critical damping, (DR)^^ = 0.7 
(p 2 + 88. 7p + 4007) A (cm) = 4007 A (in) 



(A -23) 



(3) One-tenth critical damping (DR)^ g j = 0.1 

(p 2 + 12. 66p + 4007) A (cm) = 4007 A (i (A-24) 

MAVSD System - undamped natural frequency w n ^p S j = 53-6 rad/sec 

(CT), , = 0.0993 sec. 

' '(vsm) 



(1) Critical Damping (DR)^ g ^ = 1.0 

(p 3 + 127. 4p 2 + 4030p + 40350) A, = (40350 + 4010p)A.. 

l cm ) (in) 

(A -25) 

(2) Seven-tenths critical damping (DR)^^ = 0.7 

(p 3 + 99. 4p 2 + 4030p + 40350)A^ cm j = (40350 + 4010p)A^ n ^ 

(A-26) 

(3) One-tenth critical damping (DR)^ g j = 0.1 

(p 3 + 22.83p 2 + 4030p + 40350)A (cm ^ = (40350 + 4010p)A (in) 

(A -27) 
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APPENDIX B 



REAC INSTRUMENTATION 



The instrumentation of the second order linear differential equa - 
tion of motion of the AVSD system for REAC simulation was accom- 
plished using standard techniques. For this reason, a discussion of 
the problem is not deemed necessary. The third order, linear dif- 
ferential equation for the MAVSD system required a somewhat differ- 
ent approach to effect correct REAC simulation due to the presence of 
the derivative of the input quantity. The method of instrumentation is 
developed below and a diagram of the basic REAC setup is shown in 
Fig. B-l. 

The MAVSD performance equation, as derived in Appendix A, is 
of the following general form 



desired pulse function input, instrument the following equation by 
standard REAC methods. (See Fig. B-l) 




(B-l) 



To obtain the solution of equation (B-l), where A^- n ^ represents the 





(B-2) 



Let 



A 



(in) 



y 



(B-3) 



3 2 

(p + a 2 p + a x p + a Q ) 
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I 
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Then from (B-l), 



A (cm) = b 0 y + b l p y (B ‘ 4) 

The first and second derivatives of equation (B-4) in operator form 
are: 

2 

p A ( cm ) = b o p y + b l p y (a) (B-5) 

p2 A (cm) = b o p2y + b l p3y (b) (B-5) 

From equation (B-2), 



(cm) 



b 0 A (in) 



3 2 

(p + a 2 p + a 1 p + a Q ) 



b o y 



The derivatives of (B--6) are. 



P A '(cm) = b 0 p y 

2 A I K 2 

p A'/ v = b n p y 
^ (cm) 0 ^ J 



Then, 



bi p y *.== -r- p A \ v 
1 r J (cm) 



b iP y = k p 



1 2 «. 



0 



(cm) 



Substituting in equations (B»4) and (B-5) gives: 

b l 

A (cm) “ A (cm) + bg P A (cm) 

b l 2 

pA (cm) ~ P A (cm) + bg P A (cm) 



(B-6) 



(a) (B-7) 

(b) (B-7) 



(a) (B-8) 

(b) (B-8) 



(a) (B-9) 

(b) (B-9) 
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(c) (B-9) 



2 

P 



^(cm) 



2 A . 

p A' 



(cm) 



3 p 3 A- 

b o P 



(cm) 



Multiply the equations of (B-9) by the appropriate coefficients from the 
characteristic equation of (B-l). 

b. 



a 0 A (cm) “ a O A '(cm) + a O b Q p A (cm) 

b l 2 

a l pA (cm) = a l pA '(cm) + a lI n P A (cm) 



(cm) 1 b, 



(a) (B-10) 

(b) (B-10) 



a 2 p2A (cm) = a 2 p2 A '(cm) + a 2 b^ P A '(cm) 



3 n 3 A- 

’o 



(c) (B-10) 



Note that in the basic REAC setup of Fig. B-l the quantities ( a o A '( cm )^ 
and (a^ A'^ cm j) are present. Then to determine (a Q A ( cm )) f rom 
equation (B-10)(a), it is necessary to multiply (aj) by a factor (k^) 
such that. 



a i k i 



0 b. 



or 



1 

Similarly for (B-10)(a) and (b), 

k 0 



a o b i 

■SFo 



a i b i 



a 2 b 0 



(a) (B- 11) 



(b) (B-l 1) 



k 



3 



a 2 b l 



(c) (B-ll) 



By the above procedure of instrumenting equation (B-2) and cor- 
rectly summing quantities at pertinent points, it is possible to obtain 
the displacement, velocity and acceleration responses of equation 
(B-l) to a desired pulse input to the system. The method avoids 
REAC differentiation of the pulse input and gives reproduceable time 
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function responses. 

Rectangular pulse generation was accomplished using a step in- 
put which was returned to zero magnitude at the end of the desired 
pulse duration by a relay. Displaced cosine pulse function generation 
was simply obtained by intrumenting an undamped second order sys- 
tem. Pulse amplitude was controlled by varying the initial condition 
on the second integrator. Pulse duration was established by a relay, 
and the undamped natural frequency was adjusted through proper 
choice of equation coefficients to give one cycle of the oscillation for 
the period of the pulse duration. Scale factors to meet loop gain lim- 
itations of the REAC were such as to make any minor deviation in 
pulse shape insignificant in their effect on system response. 
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APPENDIX C 

TABULATION OF NUMERICAL VALUES OF THE RESPONSES FOR 
IDEAL LINEAR AVSD AND MAVSD SYSTEMS 

Table C-l Values of controlled-member displacement ratio for rect- 
angular pulse function input to AVSD system (PD)(NP)R = 0.1 



T(NP)R 




DR = 0.1 



DR = 0.7 



DR = 1.0 



0.000 

0.050 

0.151 

0.252 

0.352 

0.453 

0.554 

0.655 

0.755 

0.856 

0.957 

1.057 

1.158 

1.260 

1.360 

1.460 

1.510 



0.000 

0.095 

0.654 

0.978 

0.915 

0.524 

0.000 

-0.472 

-0.714 

-0.671 

-0.388 

0.000 

0.363 

0.518 

0.488 

0.287 

0.157 



0.000 
0.090 
0.462 
0.523 
0.420 
0.267 
0.135 
0.045 
0.000 
- 0.021 
-0.024 
-0.019 
-0. Oil 
-0.006 
- 0.002 
0.000 
0.000 



0.000 

0.079 

0.390 

0.420 

0.333 

0.235 

0.155 

0.099 

0.061 

0.037 

0.022 

0.014 

0.009 

0.006 

0.004 

0.003 

0.002 



A 



ref 



A 



cm(max)(DR - 0. 1) 
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Ta ble C-2 Values of controlled- jnember velocity ratio for 
pulse function input to AVSD system (PDJJNPJTl = 0.1 



T(NP)R 


A /A , 

cm' ref 


DR = 0.1 


DR = 0.7 


DR = 1.0 


0.000 


0.000 


0.000 


0.000 


0.053 


0.553 


0.476 


0.430 


0.101 


0.994 


0.659 


0.566 


0.151 


0.766 


0.324 


0.213 


0.252 


0.203 


-0.068 


-0.081 


0.352 


-0.368 


-0.216 


-0. 149 


0.453 


-0.741 


-0.224 


-0.138 


0.554 


-0.803 


-0.167 


-0.103 


0.655 


-0.568 


-0.100 


-0. 072 


0.755 


-0.147 


-0.047 


-0.046 


0.856 


0.259 


-0.138 


-0.029 


0.957 


0.535 


0.003 


-0.016 


1.057 


0.589 


0.010 


-0.013 


1.158 


0.420 


0.010 


-0.008 


1.260 


0.110 


0.008 


-0.006 


1.360 


-0.184 


0.007 


-0.005 


1.460 


-0.380 


0.000 


-0.002 


1.510 


-0.430 


0.000 


-0.001 



A ref A cm(max)(DR = 0.1) 
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Table C-3 Value of controlled-member acceleration ratio for rectan- 



gular pulse function input to AVSD system (PD)(NP)R = 0. 1 



! 

T(NP)R 


A /A . 
cm' ref 


DR = 0.1 


DR = 0.7 


DR = 1.0 










0.000 


0.000 


0.000 


0.000 


0.025 


1.000 


0.810 


0.797 


0.095 


0.760 


0.334 


0.228 


0 121 


-0.334 


-0.633 


-0.659 


0.201 


-0.494 | 


-0.337 


-0.258 


0.302 


-0.494 


-0.122 


-0.006 


0.403 


-0.314 


0.000 


0.005 


0.503 


-0.051 


0.057 


0.022 


0.605 


0.208 


0.057 


0.020 


0.705 


0.357 


0.051 


0.009 


0.806 


0.359 


0.033 


0.004 


0.907 


0.235 


0.020 


0. 000 


1.007 


0.038 


0.006 


0.000 


1.108 


-0.152 


0.000 


0.000 


1.208 ~- 


-0.258 


0.000 


0.000 


1.310 


-0.268 


0.000 


0.000 


1.410 


-0.177 


0.000 


0.000 


1.510 


-0.033 


0.000 


0.000 



^ref ^cm (max)(DR = 0.1) 
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Table C-4 Values of controlled-member displacement ratio for rec 
tangular pulse function input to MAVSD system (PD)(NP)R = 0.1 



T(NP)R 




A /A , 

cm' ref 




DR = 0.1 


DR = 0.7 


DR = 1.0 


0.000 


0.000 


0.000 


0.000 


0.063 


0.152 


0.095 


0.101 


0.101 




0.259 


0.127 


0.126 


0.514 


0.368 


0.346 


0.189 


0.817 


0.535 


0.471 


0.246 


j 




0.498 


0.252 

0.258 

0.289 

0.314 

0.315 

0.327 


0.976 

1.000 

0.967 


0.583 

0.554 


0.467 


0.378 


0.837 


0.483 


0.411 


0.503 


0.270 


0.303 


0.288 


0.630 


-0.339 


0.147 


0.185 


0.693 

0.755 


-0.556 

-0.664 


0.041 


0.104 


0.794 

0.882 


-0.679 

-0.560 


-0.001 


0.050 


1.007 


-0.173 


-0.038 


0.016 


1.132 


0.263 


-0.043 


-0.007 


1.126 


0.462 


-0.041 


-0.018 


1.285 

1.385 


0.490 

0.405 


-0.034 


-0.026 


1.510 


0.118 


-0.027 


-0. 028 




~~K 7 = 

ref 


Lft 

cm (max)(DR = 


0.1) 
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Table C-5 Values of controlled-member velocity for rectangular pulse function input to MAVSD 

system (PD)(NP)R = 0.1 
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(continued on next page) 



Table C-5 - con't. 




n 



• C 
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ref cm(max)(DR = 0.1) 



Table C-6 Values of controlled-member acceleration ratio for rec- 
tangular pulse function input to MAVSD system (PD)(NP)R = 0.1 



T(NP)R 


A / A . 

cm' ref 


DR = 0.1 


DR = 0.7 


DR = 1.0 


0.000 


0.000 


0.000 


0.000 


0.013 


0.405 


0. 127 




0.025 


1.000 




0.200 


0.045 




0. 279 




0.063 




0. 206 


0.122 


0.076 




0. 110 




0.101 


0.716 


-0.022 


0.064 


0.113 




-0.178 




0.126 


-0.378 






0.138 






-0.163 


0.189 


-0.505 


-0.106 


-0.092 


0.252 


-0.550 


-0.063 


-0.047 


0.314 

tt 


-0.507 






0.378 


-0.394 


-0. 014 


-0 006 


0.503 


-0.047 


0.007 


0.004 


0.630 


0.270 


0.012 


0.007 


0.755 


0.380 


~ 


0.006 


0.882 


0.270 


0.007 




1.007 


0.027 
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Table C-ll Values of controlled- member velocity ratio for displaced cosine pulse function input 
to MAVSD syste m with (PD)(NP)R = 0.1 
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Table C-12 Values of controlled-member acceleration ratio for displaced cosine pulse function 

input to MAVSD system with (PD)(NP)R = 0.1 
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Table C-13 Values of contro lled-member displacement ratio for displaced cosine pulse function 

input to AVSD system DR = 0.7 
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Table C-14 Values of controlled-member velocity ratio for displaced cosine function input to AVSD 

system DR = 0.7 
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Table C-15 Values of controlled-member acceleration ratio for displaced 
cosine function input to AVSD system DR = 0.7 
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ref “ cm(max)[ (PD)(NP)R = 0. 1 ] 



Table C-16 Values of controlled- member displacement ratio for displaced 
cosine pulse function input to MAVSD system with DR-0. 7 
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(max) 



Table C-17 Values of controlled-member velocity ratio for displaced 
cosine pulse function input to MAVSD system with DR=0. 7 
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(ref) cm(max)[ (PD)(NP)R - 0. 5] 



Table C-18 Values of controlled-member acceleration ratio for displaced cosine 
function input to MAVSD system with DR = 0.7 
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(ref) " cm(max)[ (PD)(NP)R = 0.1] 



APPENDIX D 



TYPICAL CALCULATION OF ERROR 



The purpose of the following computation is to illustrate the ap- 
plication of the methods of Chapters 3 ard 5 to determine the average 
response error resulting from backlash. Tables 5-4 and 5-5 sum- 
marize the results of similar calculations. Figures 5-19 and 5-20 
show the variation of the average percentage error with pulse strength. 



1. Pulse Function: 



MAVSD System 
Displaced Cosine 

(PD)(NP)R = 1.0 



2 . 



3. 



Motor Limitations: At the torque summing member^maximum 

acceleration = 9300 rad/ 2 

/sec 

maximum velocity = 415. 0 rad/sec 



Gear Train: 



L (tsm) y K 



= 100 



(cm) 



Total backlash = 0. 5 milliradian = A 
Effective backlash = 0.25 milliradian 



(bl) 

' A (bl)(cm) 



4. Pulse Strength: 

(a) From Table 5-3, a pulse height of 1 mr. gives 

A, . = 1441 mr / 2 

(cm) / sec 



(b) From motor characteristics and gear ratio. 
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A, w > = 93000 rar./„ 2 
(cm)(max) / sec 



(c) (PH) 



93000 



max 1441 



- 64. 5 mr. 



(d) (PS) = 



(PH;(T n ! (64. 5)(.0993) 



3. 21 mr. sec. 



5. Peak amplitude of displacement response: 



From Table 5-3, 



A (cm)(max) < PH > x °- 967 



L (cm)(max) = 64. 5 x 0. 967 = 62. 4 mr. 



6. Peak amplitude of velocity response 
From Table 5-3, 



A, x = (PH) x 25 = 64.5 x 25 = 1612 mr./sec. 
(cm) / 



From motor characteristics and gear ratio, 



A (cm)(max) ” 4150 mr / sec - 
Therefore the velocity is well below it's limit. 

7. Application of the effective backlash to the velocity response 
of the mean point. ( Shaded triangular areas of Fig. 5-17. ) 



A (bl)(cm) 



__ AVel x Time 
2 



0. 25 mr. 



Note: As the velocity peaks decrease in amplitude, a great- 
er period of time is required to take up the effective 
backlash. 
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8, Deviation due to backlash Fig. 5-16 

(a) During time the backlash is taken up the deviation rises 
parabolically. 

(b) Deviation remains constant at the value of the effective 
backlash until the velocity reaches a negative maximum. 

(c) Deviation falls parabolically to zero as backlash is 
taken up in the opposite direction, 

9. Average error in percent 




>00 



2 



APEM 4 



x 100 




where (£)f(t.) = ordinate of deviation curve. Fig. 5-18 

f(t) = ordinate of displacement curve. Fig. 5 16 



(a) Integrating graphically. 




| f(t)|^ = 173 (mr)^ sec 




,co 



2 




10. Average error for uncertainty level (U)f(t) of 0. 5 mr. 



| f(t)||(U)f(t)|dt = 4.18 (mr) sec 



APEM 
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GLOSSARY OF SYMBOLS AND NOTATIONS 



( 5 1 

Self-defining notation developed by Draper et al' ' has been used 
throughout this thesis. 



a - 


amplitude 


A - 


angular displacement 


A - 


angular velocity 


A - 


angular acceleration 


ADFR - 


- angular displacement frequency response 


ADTTR - 


- angular displacement transient time response 


* APEM 


- average percentage error in magnitude 


” AVSD ■ 


- angular velocity signal damping 


4 ■ 


- non-dimensional frequency ratio 


bl - 


backlash 


bl(cm) 


effective controlled-member backlash 


BR ■ 


- backlash ratio 


(C) 


- correction 


cm - 


- controlled-member 
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DAR - 
DCPF - 
DRA - 
DR - 



e 



eq - 

(£)f(t) - 
F(/) - 
f(t) - 



fr - 



FR - 
I - 



id - 



IFF - 

j " 

MAVSD - 
M - 



dynamic amplitude ratio 
displaced cosine pulse function 
dynamic response angle 
damping ratio 
natural logarithm base 
equivalent 

error in the fu nction of time 

working variable form for frequency response 

working variable form for transient time response 

frequency response 

frequency ratio 

inertia 

identified 

input forcing function 
imaginary component 

modified angular velocity signal damping 
torque 
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mr 


- milliradian 


n 


- frequency in cycles per second 


n f 


- forcing frequency in cycles per second 


n ref 


- reference frequency in cycles per second 


OR 


- output response 


P 


- derivative operator 


PF 


- performance function 


PD 


- pulse duration 


(PD)(NP)R 


- pulse duration natural period ratio 


PH 

A 


- pulse height 


'* PS - 


• pulse strength 


ps - 


■ positional servomechanism or physical system 


PSR - 


• pulse strength ratio 


q - 


quantity 


REAC 


- Reeves Electronic Analogue Computer 


S - 


sensitivity - dimensional similarity factor 


t - 


time 
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undamped natural period 



T - 
n 



(T)(NP)R - 



T /T - 
p' n 



tsm - 



ttr - 



(U)f(t) - 



r - 

<p . 



time natural period ratio 
working variable expression for (PD)(NP)R 
torque summing member 
transient time response 

uncertainty in the transient time response 
angular forcing frequency 
angular undamped natural frequency 
working variable form for DAR 
working variable form for DRA 
working variable form for S 
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